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Abstract. We propose a new method to detect observational appearance 
of Dark Matter axions. The method utilizes observations of neutron stars 
(NSs) in radio. It is based on the conversion of axions to photons in strong 
magnetic fields of NSs (Primakoff effect). Whether the conversion takes place, 
the radio spectrum of the object would have a very distinctive feature - a 
narrow spike at a frequency corresponding to the rest mass of the axion. For 
example, if the coupling constant of the photon-axion interaction is M = 10 10 
GeV, the density of Dark Matter axions is p = 10 _24 gcm~ 3 , and the axion 
mass is 5 /xeV, then a flux from a strongly magnetized (10 14 G) NS at the 
distance 300 pc from the Sun is expected to be about few tenths of mJy at 
the frequency ~ 1200 MHz in the bandwidth m 3 MHz. Close-by X-ray dim 
isolated neutron stars are proposed as good candidates to look for such radio 
emission. 



1 Introduction 

Most of the matter in the Universe is "dark", i.e. it can not be observed 
directly by astronomical observations, as the particles that form it are not 
baryons [T], [2]. A lot of types of new particles were suggested by the theorists 
to explain the Dark Matter (DM) problem. One of the best candidates with a 
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strong theoretical background is axion - a light neutral pseudoscalar particle, 
which appears in spontaneous breaking of Peccei-Quinn symmetry [3] as a 
solution of the strong CP problem. 

In fact, axions are not really dark. In external electromagnetic field they 
can couple to virtual photons and produce real photons (so-called Primakoff 
effect |U) 

a + 7virt -> 7 (1) 

Several experiments that utilize the Primakoff effect are underway now. 
One direction of such studies is to look for an effect related to conversion of 
high energy axions from interiors of the Sun in strong magnetic fields in a 
laboratory [5], [6]. Also solar axions can be searched in X-ray observations 
because of conversion of these particles in the Earth magnetic field [7j,|8j. 
All these experiments are searching for "hot" and "young" axions, not "cold" 
cosmological particles originated in the Early Universe. Here we propose an 
astronomical method to detect emission due to conversion of cosmological 
axions. Direct experimental searches for such axions are in progress now in 
laboratories, too [9|. 

The strongest magnetic fields known can be found in the surroundings of 
neutron stars (NS), so under certain conditions axions that constitute cos- 
mological DM can experience the Primakoff effect in magnetospheres of NS. 
Possible use of such super-strong magnetic fields for searching light pseudo- 
scalar bosons in processes of photon conversion was previously discussed in 
several papers [TO], [Tl] . In what follows, a possibility of detection of photons 
from such conversion is studied. 



2 Theoretical model 

The axion-photon coupling is given by the following term in a lagrangian: 

Ly* = -^F^F, U 4>, (2) 

where - the axion field strength, F^, F^ u - tensor of the electromagnetic 
field strength and its dual tensor respectively, Jjj ~~ the coupling constant. 

The axion rest mass is assumed to be small due to cosmological and 
astrophysical constraints 10~ 6 eV < m a < 10~ 2 eV, the coupling constant is 
small either: M > 10 10 GeV [32], [33]. 
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The conversion probability in transverse magnetic field B is [14] (we use 
Planck system of units where energy has inverse length dimensionality in 
eqns. ©, (US}, {II])): 



where q is the axion-photon momentum difference: 




(3) 



Q = 




(4) 



m 7 is the plasma mass of a photon: 



m 7 = 0.37 -y/ n/10 8 cm~ 3 /^eV, 



(5) 



m a , i? a - the rest mass and the energy of the axion respectively. 

We estimate the probability of the conversion of DM axions to photons in 
magnetosphere of a NS using several simplifying assumptions: 1) velocity of 
a NS relative to axions is perpendicular to its rotational axis (dispersion of 
axions velocities can be neglected due to its small value [9]) 2) r-dependence 
for the magnetic field of a NS can be described by the following relationship 
B(r) = -Bo^o/?" 3 ), w h ere A) - magnetic field strength on the surface of a 
NS, r - radius of a NS. 

Flux from the axion-photon conversion changes with a period equal to a 
half of the spin period of a NS because the dipole axis is perpendicular to 
the axion flow twice during one revolution. 

If the plasma mass of photon is equal to zero the conversion is severely 
suppressed. However, density of charged particles in a NS magnetosphere is 
quite high and this makes the conversion possible. For our estimates we use 
the Goldreich- Julian density. There are claims that the plasma density in 
the region of closed field lines of highly magnetized stars - magnetars - can 
exceed the G-J density by several orders of magnitude |l5j. This is related 
to the fact that magnetars have hard tails in their spectra, discovered thanks 
to observations aboard Integral satellite ([Ej and references therein). On 
the other hand, X-ray Dim Isolated Neutrons Stars - also known as "The 
Magnificent seven", see below - which we discuss in this paper as prominent 
candidates do not have such hard tails, so we suppose that one can securely 
set the value of the plasma density in the case of the Magnificent seven equal 
to the G-J density. [17]: 
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n GJ = 7- 10- 2 -cm- 3 (6) 
T - spin period of a NS (in seconds), B - magnetic field (Gauss). 

n(r) = a\.B{r)T~ x = a^o^r^T' 1 , 
ai = 7- 10~ 2 s • cm~ 3 • G _1 

So, the photon plasma mass depends on radius as follows: 



m 7 (r) = d2\ nlr) = a 



1/2 D l/2 3/2 -3/2/71-1/2 



Here a 2 = 3.7 ■ KT 11 cm 3 / 2 eV. The conversion takes place when the photon 
plasma mass coincides with the rest mass of the axion, m a = m 1 . 

The critical radius r c and the critical magnetic field B c {r c ) can be derived 
from the following condition: 



m 7 = m a 



and therefore: 



a 



¥ 3 arB^T-^m-^r (7) 



B c = a^a 2 2 mlT (8) 
For cold axions (E & m a , [9j) the certain relation holds: 

\mf — mi\ 

q = '— |m 7 - m a | = Am (9) 



2K 



a 



The difference between the photon plasma mass and the axion mass de- 
pends on the length of axion's path L near the critical point (the conversion 
radius) before conversion takes place: 



Am 



dm y (r c ) 



dr 



3m 7 (r c ) ^. 
2rv 



On the other hand, that length L can be determined from the condition 
of maximum probability j3]): 

qL = 7T (10) 

As a result we have: 
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After rewriting the expression for the probability of conversion P 7 = 
B 2 c M- 2 q- 2 with (fH), (51). (TTOjh we obtain: 

P 7 « 20 G-WW 3 |^ 5 V 10/3 Bf T^m^M" 2 , (12) 

where coefficient 20 G~ 2 cm _1 eV 3 is needed to adjust eq. ([3]) written in Planck 
system of units to system that is used in conclusive calculations. The full flux 
of photons also depends linearly on the critical radius r c because the amount 
of axions that propagate through the region of active conversion increases 
linearly with growth of the critical radius. So, finally we can write down our 
estimate for the amount of energy that comes from the conversion in one 
second, E: 

E ~ a^ /3 a^ /3 B 2 /3 T^ 3 r 2 ml/ 3 M- 2 (13) 

It is clear that the probability rises sharply with growth of the NS spin 
period and the rest mass of axion. However, there are limitations on and 
r c (it can not be less than radius of NS ), so we should specify our candidates 
for observations for future estimates. 



3 Possibility of observation 

For subsequent estimates we assume the range of axion rest masses 0.1 /xeV < 
m a < 10 /zeV. 

Using j6j) and (JHJ) we can obtain the following relation: 

Bc = 10 10 T (^-) G (14) 

After substitution of typical values of axion rest mass and NS parameters 
into the equation for conversion probability (PT2l it can be seen that relevant 



1 Soft Gamma Repeaters (SGRs) are neutron stars with the strongest known magnetic 
fields, B ~ 10 14 " 15 G. Fields of the most magnetized pulsars do not exceed fewxlO 13 G. 
There is no evidence of existence of NS with fields stronger then 10 15 G (although the 
physical limit of field strength is about 10 18 G). 



conversion takes place only for magnetic fields B c > 10 11 G - so we need a 
NS with strong magnetic field. Also it is favorable to have a NS situated 
close to the Solar system with a spin period as large as possible. 

In our opinion, the best candidates to produce observable signal due to 
axions conversion are X-ray dim isolated NSs. Seven objects of this type are 
known (they are dubbed the Magnificent Seven (M7), and we use this term 
below), see a review on isolated NSs in p2]. They possess very strong mag- 
netic fields (up to 10 14 G) and they are located not far from the Solar system 
(~ 300 pc[19]EI). Their present evolutionary state is not known, yet: we do 
not know whether the M7 sources are similar to normal pulsars or whether 
they have already passed this evolutionary stage, and so do not produce 
radio pulsar-like emission. We use the Goldreich- Julian density of charged 
particles in magnetosphere for our estimates. If the density is significantly 
higher than the Goldreich- Julian value (about such possibility see, for exam- 
ple, [21j), the conversion would be greatly depressed, see eq. (|T2j) because it 
would take place in area with low magnetic field strength and therefore the 
probability of conversion would be tiny. 

So, let the magnetic field on the surface of NS be equal to 10 14 G, spin 
period is equal to 10 s and the distance from the Earth is equal to 300 pc. 

For m a = 5 /ieV: 

B c » 2.5 ■ 10 12 G 



3.4 r 



g 2 = ^ = 1.3-10- 16 eV 2 
2rv 



P 7 « 0.2 



It is essential to estimate total mass of axions that fly through the zone of 
active conversion (r < r c ) per unit time to obtain estimates for the energy of 
electromagnetic waves from the conversion |22| . 



m = 2vr(r c - r )GM m p V -\ (15) 



2 The closest object RX J1856.5-3754 is located at ~ 170 pc[20] 
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where p -axion density, v - NS velocity relative to DM, M NS - NS mass. 
Density of the DM, p (and therefore the density of axions) is set equal to 
10~ 24 gcm~ 3 [9] (the density may be smaller as such large values can appear 
only if a source is in a caustic, but as the flux depends on it linearly, it can be 
easily recalculated with any value of the axion density), velocity was set equal 
to v — 100 km s -1 and we use mass M NS = 1.4M and radius r = 10 km. 
Total mass of axions propagating through the "active" region per second is: 

rh = 2.8 • lO^s" 1 

The energy that comes from the conversion every second can be estimated 

as: 

E = P~,mc 2 = 5.4 • 10 22 ergs _1 

The electromagnetic flux at 300 pc from a source might be 5-10 -21 erg cm -2 s _1 . 
Radio waves would have frequencies near central frequency corresponding to 
the rest mass of axion m a , fo = 1200 MHz in the bandwidth 5f = foq/m a = 
2.8 MHz. The density of that flux might be equal to 0.2 mJy. 

We made estimates for the values of axion rest mass between 0.1 peV 
and 10 peV. Probability of the conversion rises swiftly with increase of axion 
rest mass and reach its saturation value P = 0.5 at m a « 7/ieV (figH]). The 
predicted observable flux has a sharp peak at that value of axion rest mass 
and then steeply decreases (fig. [2]). The bandwidth of the signal from the 
conversion smoothly increases with increasing axion rest mass (fig. [3)) 

The estimates above are only upper limits. We do not take into account 
the dipole structure of the magnetic field: 

„, , 3n(nm) — m . . 

£(n,r)= V J , (16) 

where n - the unit vector along the radius vector r, m the vector of 
the magnetic dipole. Accurate estimates with an exact configuration of the 
magnetic field, can slightly reduce our predictions for the flux. Also, the 
flux might be smaller because of the reverse conversion, but these questions 
require a detailed study of individual cases, and so we do not discuss them 
in this note. 

It is necessary to mention the possibility of photons absorption due to 
their propagation in high-density plasma near a NS surface. Indeed, the 
EM waves of a certain frequency can not propagate through regions where 
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plasma frequency exceeds the EM wave frequency. In the case of a NS that 
moves through an axion flow away from the Earth, converted photons will be 
reflected by the magnetosphere so the sought signal will be weakened. In the 
opposite case, dense plasma will act as a mirror, sending converted photons 
to the Earth, thus boosting the signal (it might be confined in a solid angle 
less than 4ir). 

Probability of conversion could be further decreased by the effect of vac- 
uum polarization [11]. However, the strength of magnetic field in a region of 
active conversion is well below the critical value 5 • 10 13 G. The variability 
of the signal depends on the angles \I/ and <3> between the spin and magnetic 
axis and between the spin axis and the velocity vector of a NS. 

In general, the flux variability with time is very complicated and should be 
studied numerically. The figure [4] represents the dependence of the amplitude 
of the flux variability A = (S max — S m in) / S max on the angle when the relative 
velocity of NS to axions is perpendicular to its axis of rotation. 

It is easy to see that the variability will be significant for the angles close 
to 90° (P > 0.5 if * > 80°) and usually it will not exceed 0.1. The period of 
variation might be twice shorter than the spin period of the NS. 

At the moment the most stringent upper limits on the radio emission from 
the M7 are given by [23]. At the frequency 820 MHz they are about 10 mJy. 
At lower frequencies the situation is less clear. The group from Pushchino 
[24] announced detection of pulsed emission from two NSs belonging to the 
M7. On the other hand, recent observations with GMRT (Joshi et al., poster 
at COSPAR-2008 and work in progress) do not confirm it. In the near future 
LOFAR observations [25] can be used to put better limits or to detect the 
signal due to axion-photon conversion. The signal from the conversion will be 
strongly depolarized because direction of magnetic field is different in various 
parts of active conversion region. 

4 Conclusions 

We suggest to use radio observations of close-by X-ray dim cooling isolated 
NSs to search for observational appearances of DM axions. Whether the 
axion-photon conversion in NS's magnetic field takes place, the radio spec- 
trum of the object might have a very distinct feature - a narrow spike at a 
frequency corresponding to the rest mass of axion. If the coupling constant 
of the photon-axion interaction is M = 10 10 GeV, the density of DM axions 



8 



is p = lCT 24 gcm~ 3 and a NS with B ~ 10 14 G is located at the distance 
of 300 pc from the Solar system, then the flux density of signal for axions 
with rest mass of 5 /xeV is as large as several tenths of mJy at the frequencies 
w 1200 MHz in the bandwidth » 3 MHz. 
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Figure 1: Probability of conversion vs. the axion rest mass. The probability, 
P, increases with the mass increase. There is a saturation at P = 0.5 because 
of the process of reverse axion-photon conversion that effectively suppresses 
further growth of P. Here and in the figures thereinafter, the magnetic field 
on the surface of NS is equal to 10 14 G, spin period is equal to 10 s and the 
distance from the Earth is set to 300 pc. 
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gure 2: Flux density of the signal that comes from the conversion of axions 
a certain rest mass. 
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Figure 3: Bandwidth of the possible signal from the conversion vs. the 
frequency of observation (axion rest mass). The bandwidth is defined by 
the magnitude of axion-photon momentum difference q(f) from the eqs. 

m,m),q 2 ~f 5/e 
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Figure 4: Dependence of the variability of the flux on the angle 9 between 
the axis of rotation and the magnetic dipole axis. 
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